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Abstract 

Within the ghost-free Analytic Perturbation Theory (APT), devised in the last 
decade for low energy QCD, simple approximations are proposed for 3-loop analytic 
couplings and their effective powers, in both the space-like (Euclidean) and time-like 
(Minkowskian) regions, accurate enough in the large range (1-100 GeV) of current 
physical interest. 

Effectiveness of the new Model is illustrated by the example of T(1S) decay 
where the standard analysis gives as(Mx) = 0.170±0.004 value that is inconsistent 
with the bulk of data for as- Instead, we obtain = 0.185 ± 0.005 that 

corresponds to aY°^{Mz) = 0.120 ± 0.002 that is close to the world average. 

The issue of scale uncertainty for T decay is also discussed. 


Introduction 

Theoretical expressions for measured quantities of hadron physics contain QCD running 
coupling as ■ The common formula (eq.(7) in Bethke[T] or eq.(9.5) in Particle Data Group 
(PDG) review[2]) suffers from the unphysical pole singularity that is an indispensable 
feature of all renormalization group (RG)sums of ultra-violet (UV) logs obtained from 
perturbation theory. In QGD, this issue becomes especially troublesome in the few GeV 
region. 

The hrst model of a ghost-free QED coupling a®^^(Q^;a) was devised long ago [3] 
on the basis of the Kallen-Lehmann analyticity in the complex plane. The price for 
the absence of ghost is non-analyticity of a) with respect to coupling constant 

a at a = 0 . 

The idea to use the Kallen-Lehmann imperative to get rid of unphysical singularities 
(like Landau pole) was transferred to QGD in the mid-90s[l] and named “analyticization”. 
It results in a regular in the low-Q^ Euclidean domain effective coupling function with a 
hnite IR limit aEiCp = 0). 

On this basis, another ghost-free construction for QGD effective coupling in the Min¬ 
kowskian domain aM{s) was introduced [5] via integral transformation. This result turned 
out to be equivalent to the result of 7r^-terms summation derived by similar means P, [7] 
in the early 80s (see also, 0)- 
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Later on, both the constructions were joined [9] by suitable integral transformations 
into the so-called Analytic Perturbation Theory (APT). An essential feature of the APT 
scheme is that Minkowskian and Euclidean counterparts for powers of usual QCD coupling 
{asiQ"^))^ form nonpower sets {2lfc(s)} and {Ak{Q^)} ■ For the fresh reviews of APT 

see [ini [H] • 

Meanwhile, ghost free APT expressions for effective couplings in the Euclidean 
and Minkowskian q;m(s) regions, as well as for their “effective powers” AkiQ"^) and 2lfc(s), 
are presented by simple analytic expressions only in the one-loop case (see eqs.(E])-([n])), 
which are not accurate enough for practical goals. 

Higher-loop APT expressions are more intricate involving a special Lambert function. 
A few years ago the hrst three of them Ak{Q^) , 2tfc(s); (/c = 1, 2, 3), sufficient for most 
of applications, were tabulated by Magradze and Kourashev [12] in the 3-loop case. 

Here, we propose simple model expressions for 3-loop APT couplings aE^Q"^) = 
Ai, aM^s) = 2li and for higher expansion functions Ak{Q‘^), 2lfc(s) connected by simple 
iterative relations. The Model depends on one additional parameter and is valid in the 
region above 1 GeV . 

The paper is organized as follows. Sections 1 and 2 contain a brief review of the main 
ideas, technique and some results of APT. In Section 3, we present our Model. Section 
4 contains revised analysis of T decay data with ag value extracted anew by means of 
APT and our Model. Special attention is paid to scale uncertainty. The last. Section 5, 
is devoted to the summary of the results. 

1 Outline of the Analytic Perturbation Theory 

To start this short overview, remind that the cornerstones of APT are the analyticity 
of coupling functions and compatibility with linear integral transformations. 

Here follows compendium of main dehnitions. The most elegant APT formulation is 
based on a set of spectral functions {pi(o')} dehned as 

Pk{z) = Im([a,(-z)]*^) (1) 

The hrst of them, pi = p{a) is just Kallen-Lehmann spectral density for the Euclidean 
APT coupling. Then, higher Euclidean (“analyticized /cth power of coupling in the Eu¬ 
clidean domain”) and Minkowskian (“effective kth power of coupling in the Minkowskian 
domain”) APT functions will be respectively given by 

A,{Q^) = AlaJ] = i r°° = R|a.] = i /^” . (2) 

TT io ^ + V TT cr 

These functions are related by integral transformation 

AM‘‘) = D|a,] = r°° 

Jo [s + Q^y 
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and its reverse. The differential relations connecting higher spectral fnnctions 


1 dPk o o / N 33 — 2nf 

~ p-T’j- = Po Pk+1 + Pi Pk+2 + • • • , Poi^f) = - — - , Pi 

k a ma Iztt 

induce analogous relations for expansion functions 


153 - 19n/ 
2^ 


•• ( 3 ) 


1 d2lk{s) 
k dins 


^ ^ /3n—l^fc+n('S) 5 

n>l 


1 dAkiQ^) 

k dlnQ2 


Pn-lAk+njQ'^), 

n>l 


(4) 


which can be used for iterative definitions. Numerically, these beta-coefficients (defined 
according to Bethke[T]) and their useful combination h = (3i/Pi are of the order of unitj0 

(3q{A T 1) = 0.6631 ± 0.0530 , /3i(4 ^ 1) = 0.3251 ± 0.0802 , 6(4 T 1) = 0.7392^°;°^?^ . (5) 


2 Main Results of APT 


One-loop case. In this case, the APT formulae are simple and elegant. Starting with the 

perturbative RG-improved QCD coupling = l/(/3o/), with the help of (fT|),(]2|) 

one arrives at the ghost-free effective Euclidean □ 




1/1 \ 
^ V7^ A^-QV ’ 


I = In 



( 6 ) 


and Minkowskian 

atA) 


1 

—— arccos 
PqTT 




L 



(7) 


APT couplings. Here, higher functions At , 21* can be defined via recursive relations (jl]) 
with only one term in the r.h.s. The second and third one-loop functions are 


A’(') = ^ (p 



1 g/ _ 1 _ ^21 \ 

2 (e^ - 1)3 ) ’ 


2ip\L) 

2t«(L) 


1 


1 

1 L 

^(L2+7r2)2 


( 8 ) 


In Section 4.2, we shall also need the fourth Minkowskian function, 

1 


2t«(L) = 


P! (A2 + 7r2)= 


(9) 


All APT functions obey important properties that are valid in the higher-loop case: 

^There is a misprint in numerical values of /3o(4 T 1) in paper |10j . 

^Note, we change the notation for arguments of the APT functions: I and s ^ L . 
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• Unphysical singularities are absent with no additional parameters introduced. 

• In the Euclidean and Minkowskian domains, QCD couplings aE{Q‘^) = Ai{Q ^), 
o:m{s) = 2ti and their “effective powers” AkiQ"^) 21a;(s) , are different functions 
related by integral operations A[ ] and M[ ] explicitly dehned in eqs. ([2]). 

• Higher functions, e.g., ([8]), (9) are not equal to powers of the hrst ones l|6 |) . l|7 f) . 

• Expansion of an observable in coupling powers {as{Q‘^))"' in the Euclidean or in 
(q!s(s))"' in the Minkowskian case is substituted by nonpower expansion in sets 
{Ak} , or {2lfc} respectively. The latter expansions exhibit a faster convergence, as 
compared to common power expansion. 


Two particular notes are in order: 

— The APT functions Ak{Q ‘^), 2tfc(s) essentially differ from common expansion func¬ 
tions {ds)^ in the low energy region, where they are regular with finite IR limit. In par¬ 
ticular, a£;(0) = aM(0) = l/Po- This behavior provides high stability with respect to 
variation of the renormalisation scheme |13j . 

— In the UV limit, all APT functions tend to their usual counterparts {dg)^ ■ The 
related small parameters for the measure of deviation are eg in Euclidean and sm in 
Minkowskian case 


liiqs/A")' 


( 10 ) 


Remark here that influence of APT contributions to QCD coupling and to observables 
is quite different in the Euclidean and Minkowskian regions. In the first case, APT 
correction ~ e^; with rise of momentum becomes less than 1% already at Q ~ 10 GeV. On 
the contrary, the influence of Minkowskian correction sm can be traced up to 100 GeV 
scale where the difference between two-loop (NLO) analytic and non-analytic coupling 
makes roughly 5%. Even in the 3-loop (NNLO) standard Minkowskian case, the effect of 
TT^ terms ~ afTT^/dg remains essentiao up to 10 GeV. 

Higher-loop case. The two-loop expressions are more complicated. Here, exact QCD 
coupling dg can be expressed explicitly 


in terms of the Lambert function W{z) defined as a solutioE@ of the transcendental 
equation We^ = z. This expression yields rather involved formulae for Ak and 21^, in 
terms of W_i. 

■^For more details see Sections 3.2 and 4 in [TO] . 

®One should be careful with choosing a proper branch W-i, see, e.g., [I4] . 
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In the three-loop case, one encounters more complications. Here, only for Fade ap¬ 
proximated beta-function, exact solution can be expressed explicitly [ 12 ] in terms of the 
Lambert function and one meets the same situation as with exact two-loop solutiorl^. 

The total picture becomes even more involved after taking into account the matching 
relation for adjusting kinematic regions with different values of the flavour number nj. 
The devised scheme [9], known as “global APT”, has been studied by Magradze and 
Kourashev at the two- and three-loop level. They calculated numerical tables for the first 
three functions 21^1,2,3 and ^ 1 , 2,3 at three values of = 350,400,450 MeV in the 
interval 1 GeV < a/s, Q < 100 GeV [H] [15], and 2 li ,2 , ^ 1,2 in the interval 0.1 GeV < 
\fs-iQ 3 GeV [16]. Unhappily, their numerical tables, as well as complicated analytic 
formulae, are not comfortable enough for QGD practitioners. Some other approximations 
proposed recently im are not also of wide use yet. 

Physical implications of APT. Re-examination of various processes on the basis of 
APT during the last decade has been performed in both Minkowskian and Euclidean 
regions. For a five-years old review, see paper [in- Here, we list some fresh results. 

In Ref. [18] . it has been shown that in using the standard perturbation theory for 
description of the pion electromagnetic form factor, the size of the NLO corrections is 
quite sensitive to the adopted renormalization scheme and scale setting. Replacing the 
QGD coupling and its powers by their APT images, both dependences are diminished and 
the predictions for the pion form factor turn out to be quasi-independent of scheme and 
scale settings. 

Applying appropriate generalization to the fractional powers [19] of APT coupling, 
the authors of [20] showed that the dependence of APT predictions for the pion form 
factor on the factorization scale was also diminished. 

Inclusive r-decay on the APT base has been studied by Solovtsov with co-authors. In 
paper [ 21 ], the stability of results with respect to the renormalization scale change was 
demonstrated. The QGD parameter value, A ^-=3 ~ 400 MeV, close to the world average 
was obtained in Refs. [211 [23] after due account of non-perturbative values of light quarks 
masses and summation of threshold singularities. 

One more curious application of APT appeared recently from the mass spectrum 
analysis of ground and hrst excited quarkonium states by the Milano group. There, it was 


argued [241 [25] that APT Euclidean coupling at the interval Q ~ 100 — 400 MeV 

should take values corresponding to A /=3 ~ 375MeV. 


®For the popular two-loop expression that represents (like eq.(9.5) in [2]) expansion of the iterative 
solution 



the corresponding Minkowskian counterpart is known [5]. However, numerically, it gives a rather 

crude approximation in the low-energy region. 


5 





3 Simple Model for 3-Loop APT Functions 

3.1 “One-Loop-Like” Model 

Our aim is to construct simple and accurate enough (for practical use) analytic approx¬ 
imations for two sets of functions 2tfc and .4,^ , A; = 1, 2, 3 . To reduce number of fitting 
parameters, one should better provide the applicability of the recurrent relations. 

To this goal, we suggest that one-loop APT expressions, eqs. (l6D,(lrl),(IHD, with modihed 
logarithmic arguments are used 

A"“'(0 = 4‘’W; ar"'(i) = a«(i.). ( 12 ) 

L* and being some “two-loop RG times”. 

Model functions ffT2l) are related by the “one-loop-type” recursive relations 

1 Amod 1 JAmod J 1 i Mtmod 

“ n/3o dh ~ n/3o dl ' dU ’ ~ n/3o dL, ' ^ ^ 

A simple expression for can be borrowed from |26], where a plain approximation 
for the two-loop effective log in the Euclidean region was used 

= 4 I f + T-4 = ! + MnyPT4P, 

Po [h 1 - exp(/ 2 ) J 

with b dehned in flTT|) . The structure of I 2 was inspired there by an idea of compensation 
of the hrst complex branch-cut of the Lambert function arising in the exact two-loop 
solution. This approximation was shown to combine reasonable accuracy in the low- 
energy range with the absence of singularities for ue ■ We extend this approach to higher 
functions in both the Euclidean and Minkowskian domains. To this goal, we change square 
root in “effective logs” ^ 2 ( 0 ) and hia) '■ VP + —>• \/W+~a7p with a, an adjustable 

parameter. It comes out from thorough numerical analysis that optimal value of the new 
parameter is a ~ 2 , while effective boundaries between the flavor regions have to he chosen 
on quark masses rUc = 1.3 GeV and mi, = 4.3 GeV just as in the MS scheme. 

That is, we formulate our Model as a set of equations (12) with (6) - (8) and 

L, = L 2 (a = 2) = L + folnvT^T^, h = kifl) = I + b\nVP + . (14) 

Here, L and I are dehned via common Aj^ values, like in (I6]),(I7]), for each of the havor 
region. This choice provides us with overall accuracy of a few per cent - see Table 1. 

Advantage of Model (|6 |) -(l8|). (fT^ . (flTD is that it involves only one new parameter, 
a = 2 with A|^ and Uf taking their usual values. 


6 












3.2 Accuracy of the Model vs data errors 

In Table 1 we give the maximal errors of our Model expressions (0), (ED,(El) in each Uf 
range obtained by numerical comparison with the above-mentioned Magradze tables in 
the interval of 3-loop ~ 350 — 400 MeV. 

As it follows from the Table, errors of Model for the hrst three APT functions are 
small, being of an order of 1-2 per cent for the hrst functions, of 3-5% for the second and 
of 6-10% for the third ones in the region above 1.5 GeV, i.e., in the rij = 4,5 ranges. 
However, its accuracy in the Uf = 3 region (above IGeV ) is at the level of 5-10 per cent. 

Meanwhile, relative contributions of typical LO, NLO and NNLO terms in APT non¬ 
power expansion for observables are usually something like 60-80%, 30-10%, and 10-1%, 
respectively (see Table 2 in Ref.|T0]). Due to this, the Model accuracy for many cases is 
dehned by that of the hrst model functions , provided that QGD contribution 

to an observable starts from one-loop contribution ~ a*. At the same time, for quarko- 
nium decays one meets the case with the leading contribution ~ af. There, the Model 
error is dehned by accuracy of the third Minkowskian function see Table 2. 

In this Table we compare our Model errors with some data errors in the low energy 
region. The last ones are taken from a recent Bethke’s reviews [1]. 

With due regard for data errors, we can now set some total margin of accuracy that 
our Model would satisfy. This margin may be chosen, e.g., as 1/3 of the experimental 
error bar, which is no less than 10% - see column 5 in Table 2). Then the accuracy limit, 
imposed upon the hrst APT functions could be about 3%, for the second function will 
be at least 3 times weaker than for the hrst one (say, 10%), and for the third ones, at 
least 6 times weaker (say, 20%). Analysing Table 1 according to these requirements, one 
concludes that it is not reasonable to use the Model below 0.5 GeV , whereas above it 
(or, in part of 3-havour and in the 4- and 5-havour regions) it is fully advisable. Now we 
may proceed to its practical application. 


4 T(1S) Decay Revisited 

To this goal, we take the T(IS') non-radiative decay. The main reason of this choice is 
the troubling disagreement,exceeding three standard deviations, between the coupling as 
extracted from this T decay and the world average. This can be clearly seen from Fig. 
9.2 of the Particle Data Group review [2]. 

An observable, which could be an apt illustration of the proposed Model should obey 
the following criteria: hrst, inhuence of the vr^ terms upon its value should be large enough, 
this observable should be measured up to a sufficient precision. The T(IS') non-radiative 
decay satishes these both. Indeed, on the one hand, the parameter cm dehned in flTUp is 
not small/l in the region (3-5 GeV) related to this decay: — 0-3 — 0.2 . 

^The importance of the terms in T decay was demonstrated recently. The rough estimate obtained 


7 



On the other hand, the non-radiative decay of IS'-state provides the best data precision 
(2.0% for the ratio of hadronic and leptonic widths [2]). Decays of 2S and 3S states give 
a poorer accuracy of this ratio. Besides, there is one more argument to re-examine T 
decay; the essential disagreement^ with the world average, just mentioned above. 

At the same time, quarkonium decays are not the best proving ground for the test 
of model expressions we have devised. Indeed, in the choice of the Model we made 
accent on the accuracy of hrst APT functions. Meanwhile, these decays are described by 
expressions ~ . Nevertheless, even with 7 per cent error (as it is seen in Table 2) the 

Model application to T(IS') decay will be interesting enough. 


4.1 T Widths 


The paper [7] was the predecessor of present work in revisiting T decay. 25 years ago, 
these authors showed the effect of analytical continuation to be signihcant for T decay 
width, resulting in a considerably larger A|^ (1.5 times) than the values obtained earlier 
by standard RG improved perturbation theory. Later on, the importance of tt^ terms 
was emphasized by Bjorken [8]. Nevertheless, to our knowledge, extraction of value 
was performed so far without due regard for the analytic continuation effects, with some 
exceptions mentioned below. 

A few words on other processes involving T, like its radiative decays and T production. 
All these have low data precision (typically, radiative widths are measured with the 10%— 
40% accuracy). For an extensive review of qq decay widths see [28] . 

The NLO ratio of hadronic and leptonic decay widths of the T S state was given 
in [29] (see also [30] ) 1^ 


i?x 


r (T ^ hadrons) 
r(T ^ e+e-) 


10(7r^ — 9)af (p) 
9 'Ka^{M-x) 


1 + 


otg (p) 


vr 


/So 2.78 





Here, I3q = 11 — (in PDG normalization) with n/ = 5 as MS scheme is used. 

Meanwhile, in the RG-invariant expression , scale p can appear only in the 

argument of QGD coupling as{iR) ■ To return eq.([T5D to the RG-invariant form, one could 
set fi = Mx 


i?x 


10(7r^ -9)«f(gx) A as(.ST) „ „ 

97ra^(Mx) \ vr 


5360[a3(sx) -f 2.30a^(sx)]; sy = . (16) 


Then, the issue of scale should be readdressed to the choice of sx in eq. flTOl) . 

in Section 4 of paper[Tn] gave AasiMy) — -1-0.012 that is more than 5 % correction to the NLO result. 

®This issue is absent in Bethke’s reviews [1], because there an as estimate is based only upon paper 
[27], where a reasonable Og value was extracted from T sum rules with Coulomb resummation taken into 
account. 

®This formula was reproduced in review [2| with error, but corrected in the last online PDG version. 











It was argued [30] that a proper value for scale n could be My/3 , due to the 3-gluonic 
mode of the T decay. Below, we consider the range close to My , ^/sy = 7 — 9 GeV , 
and discuss the scale uncertainty in the dual error estimate. The QED coupling a{My) 
is hxed on the T mass. 

4.2 Reevaluation of Aqcd from T Decay 

Calculations. The hrst attempt to re-evaluate extracted from T decay by proper 
taking into account analytic continuation effects was madJ^ in [^. 

Analysis analogous to [7] is performed here, with the APT expansion 

= 5360 [2l3(L) + 2.30 Ql^L)] (17) 

instead of formula ffT 6 |l which, in turn, within our Model, looks like 

= 5360 + 2.30 2lf . (18) 

By these formulae we extract A^ and ag values from fresh CLEO IIIjST] data 

RCLEO ^ 37 3 ^ g_75 

A few words on the fourth APT function . In the r.h.sides of eq. ([T7|l we use unpub¬ 

lished yet numerical results [32] and in ([T 8 ll - the Model expressions, eqs. (|H|), (Ej), (|T^ 
with logarithmic argument flTTl) : 21 ^ 4 '^(L) = 2 t 3 ^ 4 (L 2 ); L 2 = L + blny/L'^ + 27i'^ . 

As it follows from a more detailed analysi^l], at ^/s ~ 7 — 9 GeV relative error of 
2lf^°'^(L) is about 6 %. At the same time, in this interval, the ratio ~ 0.16. 

This means that the error due to the 2 nd term in the r.h.s. of eq. ffTSj) . is about 3 % and 
the total Model error for Ry°'^ about 9%. In turn, this gives 3 % Model error for Ug that 
is equal to Aas{My) = ±0.005. 

Now, numerical analysis of the ’’Exact” case yields A 5 = 210±5 MeV, that is ag{Mz) = 
0.118±0.0005 . Further on, we shall neglect by this small error. For the ’’Model” case with 
A 5 = 235, MeV, that is ag{Mz) = 0.120, one has the Model error Aag{Mz) = ±0.002. 

In line 1 of Part I, “PDG, IS”’, standard PT results on T(IS') decay are given. We 
present them not exactly as they were published in [ 2 | but recalculated along with modern 
experimental data. Line 2, marked “PDG, Fit” gives the published world average value 
described by the curve on Fig. 9.2 in [2], within the error bars of all the processes. Column 
“^(Mt)” means ‘As, calculated at the mass of T, according to eq.(9.5) of [2]”. 

^'^One should remark here that some other kind of summation, namely of l/u—type terms, for the im¬ 
provement of extraction from T production cross-sections was devised in |27j . Unlike the present case, 

where we deal with Ry for T decay, the authors of m studied sum rules for ratio R{s) = ■ 

Analytic continuation effects were not taken into account. 

^^Based upon calculation [32| of exact APT formula. 
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Line 1 of Part II, “Exact APT”, presents resnits of TIS decay caicniated by exact 
nnmeric tables for 2 I 3 and by [32] for 2 I 4 APT fnnction. Line 2, [Mod], presents valnes 
obtained from T(IS') decay data by means of the Model eqs.(12),(14). Here, model errors 
combine Model errors of both the terms in the r.h.s. of eq.(18). Line 3, “Crnde APT” 
gives an earlier resnlt [TO] with approximate APT estimate used to correct (Corr) the 
Bethke-2000 value as{My) = 0.170 extracted there from all the T decays data. 

Everywhere in Part II, we translate our results from APT expressions, by use of related 
A 5 values into standard non-analytic coupling at Mx and Mz for ease of comparison 
with Part I and other standard sources. 

The importance of taking into account the analytic effects for proper as extraction in 
the low-energy range is clearly seen from the Table 3. It is also evident that the Model 
error is small enough even in the non-favorable case of quarkonium decay, in particular 
when compared with the scale error. 

The scale uncertainty. Besides “direct” data error that in our case is 2% for i?x , there 
is a “hidden” uncertainty in theoretical equations flTTI) and flTSD related with a choice of 
argument L = ln(sx/A^). In given hgures, we use sx = . To discuss effect of the sx 

variation, return to issue mentioned above at the end of Sect.4.1. . 

This reference scale issue is actual for QCD analysis of all low energy data. The RG 
non-invariant term ln(Mx//i) in the r.h.s. of eq. fllSI) just represents an attempt to take 
into account this effect. The scale effect enlarges quickly with the uncertainty rise. E.g., 
for its value 1 GeV, i.e., for 8-9 GeV interval, we have 205 < A 5 < 240 MeV, while to 
the 7-9 GeV case — 180 < A 5 < 240MeV. In the Table 3, where the APT results are 
compared with ones of standard PT, for the scale error [...j^c we conditionally give hgures 
related tc@ the 2 GeV interval. 

Generally, the scale issue is an intrinsic problem of renormalization group application 
to observables. On the one hand, the common, “vulgaris”, version of RG algorithm 
corresponds to the UV (massless) case with simplihed dehnition of effective (running) 
coupling a{q^/jj^). Here, it is tacitly assumed that vertex function V{q\,ql,ql) entering 
into the a dehnition is taken with equal arguments. On the other hand, the condition 
is used. Both the assumptions are not valid for the low energy QGD case. On 
these items, the reader could be addressed to mass-dependent RG formalism [33[ [3T[ l35]. 

5 Conclusion 

1 . Our theoretical result is model explicit expressions flT^ for the analyticized 3- 
loop couplings ([6|),([7]) and their ehective powers ([8|),([9|) in both the Euclidean and 
Minkowskian regions. These are just one-loop analytic expressions of Analytic Per¬ 
turbation Theory with modihed logarithmic arguments ffT21) . The accuracy of the 

the case Vsx = My/2 = Mb there corresponds ae{Mz) ~ 0.114. and to Vsx = My/S [301 “ 

~ 0 . 110 . 
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Model is estimated to be sufficient for practical purposes in the region 1-100 GeV, 
which hosts many important processes. 


2. To illustrate the APT and Model application, we consider the as value extraction 
from data for T(1S) non-radiative decay measured with 2% error. Fixing the scale 
at the T mass, we got (for the coupling) by exact APT numerical calculation 
and by our APT Model 

«rTe""ct(Mz) = 0.118(l)e.p , af™°^^\Mz) = 0.120{2)Mod. (20) 

The comparison with result of T(1S) usual analysis, as{Mz) = 0.112(2), and with 
the world average, as{Mz) = 0.1185(20), conhrms the validity of the devised Model. 

3. It was established that the scale uncertainty essentaially reduces the value of theo¬ 
retical analysis of T decay. In the considered case, the scale error dominates over 
experimental and Model ones. This issue is worth urgent further examining. 
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Table 1: Maximal errors of the Model with a = 2 (old) 



ErOii 

Er^2 

i^rSls 

ErAi 

ErA2 

ErAs 

3 (0.5 - l.OGeV) 

3% 

5 

15 

3 

10 

20 

3 (1 - 1.5GeV) 

3% 

4 

8 

2 

7 

10 

4 

1 % 

3 

8 

2 

10 

20 

5 

1 % 

3 

6 

0.5 

4 

15 


Table 2; Maximal errors of the Model with a = 2 (corrected) 



Er^i 

Er^2 

Er^3 

ErAi 

ErA2 

ErAs 

3 (0.5 - l.OGeV) 

3% 

17 

39 

5 

16 

8 

3 (1 - 1.5GeV) 

4% 

5 

9 

3 

7 

7 

4 

0.5% 

3 

9 

1 

10 

21 

5 

1 % 

3 

6 

0.5 

4 

14 


Notation: chan. - t = Euclidean or s = Minkowskian channel; scale q = ^/s or q = \f ^; 
^aja. - relative error of a at a given scale; DIS (Bjork.) - Bjorken scaling violation in DIS, 
GLS — Gross - Llewellyn-Smith sum rule, and * marks combined theor. and exper. errors. 
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Table 3: Results of various Ug extraction from Upsilon decays 


Part I. Non-APT treatment 

Source 

cr(iVfx) 

a{Mz) 

A!V_=^ 

MS 

PDG, T,l^ 

0.170(4) 

0 .112(2) 

1461]? 

PDG, global Fit 

0.182(5) 

0.1185(20) 

217li 

Part II. APT treatment 

Exact APT, IS 

0.1805(12)e.p 

0.1179(5)e.p 

210(5) 

[Mod], 15 

0.185(5)m 

0.120(2)m 

235(25) 

[To] Crude APT 

0.183 

0.119 

222 

Exact ± [scale], 15 

0.180[7],e 

0.118[3],e 

210 [30] 
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